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The differential equations describing the flap-lag-torsional motion of a 
flexible rotor blade including third-order non-linearities were derived in [1] 
for both cases of nover and of forward flight [equations (5.2a-d)]. Making use 
of the two boundary conditions G^(x = l, t*t) = 0 and u(x»=0, t-t)*0, 
those equations were reduced to a set of three integro partial differential 
equations written in terms of the flexural deflections v(x,x), w(x,t) and the 
torsional variabl “(x,"^). These latter equations are of the form (equations 
(5.10) to (5.12) in [1]) 

G' = h (a = v,w) (1) 

a a 

= Q,-. 

X X 


For a blade with the end x = 1 free, G^ (x = 1, i = t ) = 0. The expressions for 

G , h (a = v,w), A„ and the aerodynamic term. Q, are given in [1]. 

a a 9 ' f 

X x 

Since the publication of [1] the case of hover is being addressed. Because 
of two of the boundary conditions for equations (1) are G (x = 1, t = t) = 0 
for a = V and a = w, it is convenient to first integrate equations (1) from x = 1 
to X = X to obtain 

X 

G = ^ h dx (a = v,w) (3) 

Cl ^ O 

1 ^ 

It should be noted that the expressions for h in equations (3) contain the 
generalized aerodynamic forces and Explicit forms for the third order 

expressions for Q and Q in terms of the deflections v and w, and their 
derivatives, were developed by making use of the implicit expressions given in 




Chapter VI of [11. 



At approximate solution to equations (2) and (3) was vrlttan in tarns of 


a set of orthogonal functions (x) and gj(x) (j ■ 1,2,...) as 

n n 

v(x,t) - I f.(x)v (t) ; w(x,t) - J f.(x)w .(t) 


j-1 




j' ' tj 


(4) 


j-1 


- I g.(x)e (t) 
J-1 J ^ 


(5) 


and equations (2) and (3) were then reduced to a set of ordinary differential 
equations for and (J ■ l,2,...,n). For this, ve found It convenient 

in terms of obtaining somewhat simpler "Galerkln integrals" and also in terms 
of comparing the present work with that of other authors such as [2], to pre- 
multiply equations (3) by f^(x) (1 - 1,2, ...,n) to obtain 


1 1 X 

/ f'(x)G dx - / f'(x) / h dxdx 

^ * a - 


(6a) 


1 lx 

/ f'(x)G dx ■ / f*(x) / h dxdx 
« 1 w 11 ; w 


(6b) 


As an Illustration of the above process consider the right hand side of 
equation (6a) where, as obtained from [1], 


e e 

h - V - V 
V 


- 2w sB - 2(c8) / tv'v' + w'w' + j v dxjdx - Q + O(c^) 
0 “ 1 


(7) 


Making use of equations (4) and (7), equation (6a) yields, for 1 ■ l,2,...,n 


/ fj(x)c^dx .8) - 


4c B 
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-! I Q^dxdx + O^r 
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where 


ij 


0 


lx 11 

f f' i f.dx = f.(01 ! f.dx + r f.f dx 
i 1 1 n J n 1 .1 


(9a) 
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M_ .. . = - I 
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Tc 0 (l*’) equation 


f! 
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X X 


1 0 


f . dx dx dx dx 
J 


(8) is of the same form as equation (B3a) in 


(9c) 


A complete set of shif ted Lepcnlr^ polynomials , defined in the interval 

1 

0 X ■■ 1 and normalised, for convenience, so that ■ f.f. dx =0 for i ^ 1 
- - 5 1 J 

and I for i = j, was chosen as the set ^f/'. In addition, normalized non-rotatinp 

torsional beam modes was taken as the set i p . 1 . The "Calerkin coefficients" that 

' 1 

appear in the reduced set of ordinary differential equations for ‘'•'^^(t), w^^(t) 

and f . ( ) were evaluated and stored on tape for use in the subsequent analysis 
1 1 

of the motion. The number of different coefficients in the equations is consider- 
ably reduced when orthogonal polynomials is used instead of non-rotating flexural 
beam nodes. Analytical expressions for the coefficients, that were evaluated 
recursively, were developed. As a result, the numerical value of each coefficient 
obtained is essentially "exact." For those coefficients that were expressed as 
integrals of a product of a polynomial in x and sin a x (a = (i - -^)"), the 
analytical expressions developed contained the sum c ' two series, one being the 
product of a polynomial in x and sin a x, and the other the product of another 
polynomial in x and cos a x . Each integration process was then transformed into 
an algebraic manipulation of these scries. 



It is worth noting that a nuiiiber of numerical difficulties affecting the 

accuracy of several of the coefficients was initially encountered when attempting 

to evaluate them numerica]ly. At first the numerical evaluation was done by 

making use of the IBM SSP subroutine DQATR. A number of coefficients were 

evaluated analyticallv in order to ascertain the validity of the numerical 

results obtained. Use of DQATR for this purpose was eventually abandoned when 

several of the numerical computations presented accuracy problems that yielded 

unacceptable results. At that point we have decided to change the numerical 

1 

scheme for evaluating anv coefficient, sav a = / g(x) dx , as a >= y(l) where 

^ ' 0 

v(x') = ( gCx) dx. The coefficient > exemplified here was then evaluated by 

0 

numerically integrating the differential equation > ' (x) = g(x), with '>(0) * 0, 
iiot:. x = 0 to X = 1. The numerical integrations were performed by using the 
IMSL routine DASCRU . This process yielded accurate results for a considerable 
number of the coefficients with which numerical problems were encountered with 
DQATR, and even required less computation time. However, it was eventually 
abandoned after we have found several coefficients determined numerically to be 
in error when a number of spot checks were made, and after several attempts to 
uniformly correct the problem by re-adjusting the integration step size failed. 
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